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Abstract. In this paper, we revisit the Cardassian model in which the radiation energy component is 
included. It is important for early epoch when the radiation cannot be neglected because the equation 
of state (EoS) of the effective dark energy becomes time variable. Therefore, it is not equivalent to the 
quintessence model with a constant EoS anymore. This situation was almost overlooked in the literature. 
By using the recent released Union2 557 of type la supernovae (SN la), the baryon acoustic oscillation 
(BAO) from Sloan Digital Sky Survey and the WiggleZ data points, the full information of cosmic mi- 
crowave background (CMB) measurement given by the seven-year Wilkinson Microwave Anisotropy Probe 
observation, we constrain the Cardassian model via the Markov Chain Monte Carlo (MCMC) method. 
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A tight constraint is obtained: 



-0.0479 



+0.0730+0.142 
-0.0732-0.148 



in 1,2<T regions. The deviation of Cardassian 



model from quintessence model is shown in CMB anisotropic power spectra at high Vs parts due to the 
evolution of EoS. But it is about the order of 0.1/% which cannot be discriminated by current data sets. 
The Cardassian model is consistent with current cosmic observational data sets. 
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1 Introduction 

Since the discovery of current accelerated expansion of our 
Universe [IJ[2] , a flood of models have been designed to ex- 
plain this late time accelerated expansion phase. For the 
reviews, please see [3^4,5,6,7,8,9 . In the phenomenologi- 
cal perspective, the accelerated expansion of our Universe 
can be realized through modifying the form of Friedmann 
equation via the introduction of an extra exotic energy 
component, dubbed dark energy, which has negative pres- 
sure, or by some possible modifications of gravity theory, 
say f(R) and brane models etc. As a result, the conven- 
tional Friedmann equation can be modified into the form 
of 

f(p), (1) 



H 2 



where f{p) is a function of energy density p which may 
include dark matter and extra energy components, and 
H is the Hubble parameter. The Cardassian model firstly 
proposed by Freese and Lewis [10] , where the Friedmann 
equation was modified into the form of 



i.e, cold dark matter, baryons and radiation. For the ori- 
gin of p n term, one can find several explanations [TT]. It 
can mimic the brane model which include a power law 
term due to the embedding of our universe into a five di- 
mensional bulk. And, the late time accelerated expansion 
of our Universe is because of the leakage of the gravity 
force at the large scale in the brane world. The Cardas- 
sian model can reduce to /1CDM model when n = and 
to conventional CDM model with vanishing cosmological 
constant when A = 0. The new term p n will dominate the 
energy component at the late epoch to provide an accel- 
erated expansion, then the values of n should be < 2/3. 



This model has been confronted by cosmic observa- 
tions extensively, when the energy density is composed of 
cold dark matter and baryon, for the recent result please 
see H2] and references therein, where n = w+1 = —0.039 



H 



8ttG 



P + A P n 



was obtained by using SN la, BAO, CMB shift parameters 
and observational Hubble data. In this case, the effective 
equation of state (EoS) of effective dark energy (the sec- 



+0.135 
-0.153 



(2) ond term of Eq. ©) is w 



eff _ 

de 



1. So it corresponds 



to explain the current accelerated expansion of our Uni- 
verse. Here p can be composed of conventional matter, 
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to quintessence dark energy model where the matter and 
dark energy are included only. However, it is not always 
true. When the radiation energy component is added in 
the Cardassian model, i.e. p = p m + p r , it is not equiva- 
lent to the quintessence model anymore due to the time 
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variable effective EoS [13] 



■w 



eff 



p r0 a~ 



3 p m oa 3 + pro®' 



(3) 



where po and pmO = Pco + Pbo are the present energy den- 
sity for radiation and matter and we have normalized the 
scale factor to oq = 1. At the early epoch after the recom- 
bination and before the last scattering surface, the radi- 
ation and the matter energy components are important, 
where the effective EoS w e JJ is not a constant. Therefore 
it is important to take the general form of EoS ([3]) into ac- 
count when one constrains Cardassian model from CMB 
observations. However, in the literature for example |12j . 
this important thing was not considered at all, with ex- 
ception of [13] (see also [2]), where the locations of the 
peaks of the CMB anisotropic power spectrum were used 
as cosmic constraint. Though the time variable effective 
EoS was considered in |13j , only the positions of peaks of 
CMB power spectra were used as cosmic constraints. In 
that paper [T3], by fixing the values of n s and giving a 
number of values of n, the authors gave the ranges of J? m o 
and h. In fact, it is not enough to give a tight constraint 
to model parameter space. In Ref. [Tl], the authors inves- 
tigated the CMB TT and matter power spectra. But, they 
did not give the model parameter space from the cosmic 
observations. 

In the literature, the Cardassian model was constrained 
by the so-called CMB shift parameters, i.e. R, l a and z* 
which are obtained based on a ylCDM model, from high 
redshifts. The problem that was almost overlooked is that 
the values of the CMB shift parameters depend on the cos- 
mological models, here ylCDM model. The potential logic 
is that since /1CDM model is a concordance model then 
any model which deviates slightly from /1CDM model is 
a competitive model too. However it is non-proper when 
one uses cosmic observational data points to constrain any 
other cosmological models. Strictly speaking, it is just a 
test of the possible viability of a model Q not a constraint 
to model. This is because that the potential circular prob- 
lem is committed and that is what we try to avoid in 
the cosmological constraint issue. In fact, the CMB shift 
parameters would be different for different cosmological 
models due to different physics process around the last 
scattering surface, for example early dark energy model 
[15] where the contribution from dark energy may not 
be neglected due to nontrivial equation of state (EoS) of 
dark energy. In this situation, it is dangerous to use data 
points derived in ylCDM model due to much departure 
from ylCDM model. One the other hand, the full CMB 
data points contain more information than the shift pa- 
rameters. For instance, at late epoch when the dark energy 
is dominated, the gravitational potential becomes evolu- 
tion. It affects to the anisotropy power spectra of CMB 
at large scale (low I parts) due to the so-called integrated 
Sachs- Wolfe (ISW) effect which is sensitive to the proper- 
ties of dark energy. Apparently, the CMB shift parameters 
do not include this information. So, one can expect a tight 



constraint to cosmological models when the full informa- 
tion for CMB is included. So, in this paper, we should take 
the observational data points from CMB directly not the 
derived model dependent CMB shift parameters. 

To fill out the gap and to avoid the so-called circu- 
lar problem, in this paper we shall use the type la su- 
pernova (SN la), baryon acoustic oscillations (BAO) and 
full WMAP-7yr data points to constrain the Cardassian 
model. 

This paper is structured as follows. In section [2] we 
give a brief review of Cardassian model and present the 
comic observational data sets and constraint methodalogy 
used in this paper. Section [3] is the conclusion. 



2 Constraint Methodalogy and Results 

At first, we give a grief review of the Cardassian model. 
After the definition of an effective dark energy p e JJ ', the 
Friedmann equation ^ can be recast into 

8ttG 



(p + pT) 



(4) 



where p = p m oa 3 + PrO a 4 and p e JJ is given by 



Pd, 



eff 



8irG F 

3-^0 (1 — ^m0 — ^ro) 



or equivalently 



eff _ eff ( ^m0 ffl 3 + ^r0 a 
de P deO I Ci \ Ci 



4 \ n 



(5) 



(6) 



where Q m o = ^co + ^bo and f2 c o = 8itGp c o/3Hq J?bo = 
8irGpbo/3H 2 and f2 r Q — 8-7rGp r o/3H q are present dimen- 
sionless energy density for cold dark matter, baryon and 
radiation respectively. Here we have defined the dimen- 
sionless energy density for effective dark energy as fi^el — 
8%Gp e JJ > /3Hg = 1 - Q mQ - Q r0 for spatially flat FRW 
Universe. One should notice that n is the only model pa- 
rameter. 

We consider the perturbation equations for effective 
dark energy in a spatially flat FRW Universe. We treat 
the effective dark energy as a perfect fluid with EoS ([3]) • In 
the synchronous gauge, using the conservation of energy- 
momentum tensor = 0, one has the perturbation 
equations of density contrast and velocity divergence for 
effective dark energy 

tide = -(1 + Wde)(0de + ^) ~ 3H(c^ - W de )5 de (7) 



Tde 



-H(l - 3c 2 s )0 de + -^k 2 S de - k 2 a de (8) 
1 + w 



In the sense of small deviation from ylCDM model. 



following the notation of Ma and Bertschinger [TB] . For the 
perturbation theory in gauge ready formalism, please see 
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[17]. The shear perturbation Ude = is assumed and the 
adiabatic initial conditions are adopted in our calculation. 

In our analysis, we perform a global fitting to deter- 
mine the cosmological parameters using the Markov Chain 
Monte Carlo (MCMC) method. We modified the publicly 
available cosmoMC package [TH] to include effective dark 
energy in the CAMB [19] code which is used to calculate 
the theoretical CMB power spectrum. The following 7- 
dimensional parameter space is adopted 



P = {w b ,w c ,0 s ,T,n,n s ,log[lO iu A s ]} 



(9) 



where uj b = f2i,h 2 and uj c = f2 c h 2 are the physical density 
of baryon and cold dark matter respectively, &s (mul- 
tiplied by 100) is the ration of the sound horizon and 
angular diameter distance, r is the optical depth, n is 
the newly added model parameter related to Cardassian 
model, n s is scalar spectral index, A s is the amplitude 
of the initial power spectrum. The pivot scale of the ini- 
tial scalar power spectrum k S Q = 0.05Mpc _1 is used in 
this paper. The following flat priors to model parameters 
are adopted: co b G [0.005,0.1], u c G [0.01,0.99], O s G 
[0.5,10], t G [0.01,0.8], n G [-1,2/3], n s G [0.5,1.5], 
log[10 1( M s ] G [2.7,4]. Furthermore, the hard coded prior 
on the comic age lOGyr < to < 20Gyr is also imposed. 
Also, the physical baryon density oj b = 0.022 ± 0.002 [20] 
from big bang nucleosynthesis and new Hubble constant 
H = 74.2 ± S.Gkms^Mpc -1 [21] are adopted. 

To get the distribution of parameters, we calculate the 
total likelihood C oc e~ x / 2 , where % 2 is given as 



2 _ 2 i 2 i 2 

X — Xcmb + Xbao ~t Xsn- 



(10) 



The 557 Union2 data [22] with systematic errors and BAO 
[23] are used to constrain the background evolution, for 
the detailed description please see Refs. [23]. SN la is used 
as standard candle. And BAO is used as standard ruler. 
To use the BAO information, we obtain the baryon drag 
epoch redshift Zd numerically from the following integra- 
tion [261 



obtained by using the fitting formula [27 for the fidu- 
cial cosmology. Here Dy(z) = [(1 + z) 2 D 2 A cz/ 'H(z)] 1 ^ 3 is 
the 'volume distance' with the angular diameter distance 
Da- In this paper, for BAO information, the SDSS data 
points from [23] and the WiggleZ data points [25] are used. 
For CMB data set, the temperature power spectrum from 
WMAP 7-year data [25] are employed to constrain the 
model parameters related to initial conditions. 

The constrained results are shown in Tab.[TJand Fig.Q] 
As a comparison to the result obtained in |12) where n = 
was obtained, a tight constraint was achieved 



-0.039 



+0.135 
-0.153 



in this paper apparently. To sec the difference between the 
Cardassian and the quintessence model, we also redid the 
same process by using the same data sets combination, the 
obtained parameter space with 1, 2er regions is also listed 
in the last column of Tab. [TJ 



Prameters 
f2 b h 2 



n(w) 
n s 

log[10 10 A] 
Age/Gyr 

Ho 



Cardassian 



n noot;+o.ouoi)2i+o.ooi04 

U.U__O_ .000525-0 .00102 
n 11/1+0.00449+0.00877 
u - 11 ^-0. 00453-0. 00872 
1 n o Q i +0.00255+0.00502 
1 - UJ91 -0. 00256-0. 00507 

n nsfis+0' 0(,(i3(i +0' 0242 

U.UOOO_o. 00698-0. 0233 
_fl A/I7Q+0. 0730+0. 142 
U.Wil »_ .0732-0.148 
fl Qfi«+0. 0127+0. 0256 
U.»UU_ .0128-0.0252 
o no-re+0. 0338+0. 0681 
O.UO ( O_0.0338-0.0G52 
fl 791 +0 0169+0.0320 
u - ' ol -0. 0168-0. 0351 
1 9 7/IQ+0. 109+0. 211 
ID. /^y_ .109-0.214 
fl ocq+0. 0168+0. 0351 
U.ZD»_ .0169-0.0320 
10 CO/I +1-204+2.346 
1U.OZ1_j .206-2.400 
71 1 nq+1. 863+3. 727 
I l.iaO_ L 864-3. 607 



Quintessence 



n noO /l+0.00051Y+0. 00103 
U.UZZ1_ .000520-0.00102 

114+0.00444+0.00885 
U.ll*_ . 00441-0. 00859 

1 n9Qn+ '0° 26 o+0'00 518 

l.UO»U_o 00262-0. 00520 
nSRS + (, '00 ( '25 + ().()25(] 
U.UODO_o. 00719-0. 0230 

_1 n_Q9+ '0 745 + ' 142 
± .U1»__ 0751 _ .i55 

n QfiK+°-o l26 + - 0254 

U.»UU_ .0126-0.0248 
o nS 77+0. 0335+0. 0l,!)l 
•3.UO! 1 -0.0337-0. 0657 

n 7Qn+ '0 167 + ' 0311 

U. ( OU_ .0168-0 .0344 
1 9 7CC+0. 108+0. 219 
lO. lOU_ .109-0.219 
n O7A+0. 0168+0. 0344 
U._ (U_ . 0167-0. 0311 
10 CQQ + 1 1869+2.447 
iu ' JOO -l .186-2 .390 
71 i m+1. 860+3. 745 
' 1 - llu -l. 877-3. 667 



Table 1. The mean values of model parameters for Cardassian 
and quintessence models with la and 2a errors, where WMAP 
7-year, SN Union2 and BAO data sets are used. 



rvo 

T (Vd) = / dri'f d 
Jn 

rZd drjx e (z)a T . 

az— = 1 

lo da R 



(11) 



where R = 3p(,/4/> 7 , ctt is the Thomson cross-section and 
x e (z) is the fraction of free electrons. Then the sound hori- 
zon is 

rv(zd) 

r s (z d )= / dr,c s (l + z). (12) 



where c s = 1/^/3(1 + R) is the sound speed. Also, to ob- 
tain unbiased parameter and error estimates, we use the 
substitution [25] 



, r s (z d ) . . 
d z ; -. , r s (z d ), 
r s (Zd) 



(13) 



©1 1^ 



where d z — r s {z,d)/ Dy(z), f s is evaluated for the fiducial 
cosmology of Ref. [23], and Zd is redshift of drag epoch 



Fig. 1. The 2-D contours with la, 2a regions 
marginalized distribution for Cardassian model. 



and 1-D 
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To understand the effect of model parameter n to CMB 
power spectra, we plot the Fig.[2]by taking different values 
of n where the other relevant parameter are fixed to their 
mean values as listed in the second column of Tab. [TJ 
Large values of model parameter n increase and depress 
the CMB power spectra in the left and right sides of the 
first peak respectively. It is on the contrary for small values 
of n. For the fixed current ratio i?^{g /i? m o, a larger value 
of n makes the effective dark energy dominate earlier, then 
the gravitational potential decays more, which enhance 
the lower Vs part of CMB power spectra due to Integrated 
Sachs- Wolfe (ISW) effect. One can also read off that large 
values of n shift the position of peaks to the low Vs part 
slightly as shown in the bottom panel in Fig. [2j which 
comes from the change of the ratio of sound horizons at 
current and last scattering surface due to different values 
of EoS of effective dark energy. 

Correspondingly, we show the comparison to ylCDM 
model in Fig. [3] where the mean values are taken from 
the second column of Tab. QJ From Fig. [3J one can find 
the Cardassian model match vlCDM model very well, and 
it almost in the center of error bars of WMAP-7yr data 
points. It means that current cosmic observational data 
points cannot discriminate Cardassian model from /1CDM 
model. 
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Fig. 2. The effect of model parameter n to CMB power spec- 
tra, where the solid black lines for the mean values taken from 
Tab. [TJ The dashed and dotted lines are for the cases of differ- 
ent values of n where the other relevant parameters are fixed 
to their mean values as listed in Tab. [TJ The bottom panel 
shows the corresponding ratios to the mean value case, i.e. 

V^l meanj/^l mean- 
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Fig. 3. The CMB Cf T 

power spectrum v.s. multiple moment 
I, where the grey dots with error bars denote the observed 
data with their corresponding uncertainties from WMAP 7- 
year results, the red dashed lines are for the Cardassian model 
with mean values as shown in Table [JJ the blue solid lines 
are for /1CDM model with mean values taken from [28] with 
WMAP+BAO+ifo constraint results. The bottom panel shows 
the ratios to ylCDM model. 



Now, we are in the position to see the effect to CMB 
power spectra due to the evolution of effective EoS at early 
epoch. At first, the evolutions of EoS of effective dark en- 
ergy are shown in Fig. 0] where different values of model 
parameter n are adopted and the other relevant param- 
eters are fixed to their mean values as listed in Tab. [TJ 
One can see the differences of the evolutions of EoS with 
respect to scale factor a at early epoch apparently. And 
we expect to see the effects on CMB power spectra. 

To do that we plotted the CMB power spectra for 
Cadassian and quintessence model in Fig. [5] where the 
model parameters are fixed to their corresponding mean 
values as listed in Tab. [TJ From Fig. [5j one can easily 
see the differences between Carddssian and quintessence 
model that appear at high I's part, i.e. the epoch around 
the last scattering surface as expected. However, the devo- 
tion is about the order 0.1/% which almost cannot be dis- 
criminated by current data sets. But on theoretical level, 
they are really not equivalent. 
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Fig. 4. The evolutions of EoS of effective dark energy, where 
different values of model parameter n are adopted and the 
other relevant parameters are fixed to their mean values as 
listed in Tab. [T] 



involved in a consistent way. We find out that the Cardas- 
sian model is consistent with current cosmic observational 
data sets. The constrained results are shown in Tab. Q] 
The results show that the model parameter n in 1, 2a re- 
gions are n = -0.0479t°;^^t°'i48 which are tighter than 
that, n = -O.OSQto;^, obtained recently in [12]. In Ref. 
|T3"] . the authors had used the positions of CMB peaks 
to constrain the Cardassian model. They found out that 
the values of n depends on the values of n s . They also 
pointed out that for n s — 1 the n = (ylCDM model) is 
not included in the allowed region. In this paper, we have 
used the full information from WMAP-7yr results not the 
0.0010 positions of peaks of CMB or the derived CMB shift pa- 
rameters in ylCDM model to constrain the model param- 
eter space. The results show that n — (ylCDM) case 
is included in the allowed region. It means that currently 
available data points from SN, BAO, HST and CMB can- 
not distinguish a Cardassian model from a ylCDM model. 
We also show the effect of time variable EoS of effective 
dark energy to the CMB power spectra for the Cardassian 
model by comparison to that for the quintessence model. 
We find out the deviation is about the order of 0(1O -3 ) 
which cannot be discriminated by current CMB data sets. 
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I 

Fig. 5. The CMB power spectra for quintessence model and 
Cardassian model, where the values of model parameters are 
adopted as listed in Tab. [1] The bottom panel shows the devi- 
ation from quintessence model. 

3 Conclusion 

In summary, in this paper we have revisited the Car- 
dassian model where the radiation energy component is 
included. In this situation, the EoS of effective dark en- 
ergy for the Cardassian model is not a constant but time 
variable. Then the Cardassian model is not equivalent to 
quintessence model anymore. We performed a global fit- 
ting on the cosmological parameters in Cardassian model 
by using a completely consistent analysis where the full 
information from WMAP-7yr released data points were 
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